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The Hubbard on-site repulsion U between opposite spin electrons on the same atomic orbital is 
widely regarded to be the most important source of electronic correlation in solids. Here we extend 
the Hubbard model to account for the fact that the experimentally measured atomic U is different 
from the one obtained by calculation of the atomic Coulomb integral. The resulting model describes 
quasiparticles that become increasingly dressed as the number of electrons in the band increases. 
Superconductivity can resuit in this model through quasiparticle undressing. Various signatures of 
this physics in spectroscopies in the normal and superconducting states are discussed. A novel effect 
in the normal state is predicted to be electroluminescence at the sample-positive counterelectrode 
boundary. 
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The Coulomb repulsion integral for two electrons in a 
hydrogenic Is orbital for a nucleus of charge Z is 
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(«o =Bohr radius). However, the effective atomic on- 
site repulsion Ueff{Z) is obtained by considering the dif- 
ference in energy between ions with different number of 
electrons. For Hydrogen (Z = 1) and Helium (Z = 2) 
the experimental values are [|l| 

Ueffil) = I - l2MeV = U{1) - A.lbeV (2a) 



C/e//(2) = /// - // = 29.92ey = U{2) - A.lQeV (2b) 

Here, / = 13.6eT/, A = O.lAleV, h = 2AA%eV, In = 
bAAOeV are the ionization energy and electron affinity 
of H , and the first and second ionization energies of He. 
Remarkably, the difference between Uef / and U is nearly 
the same ioi Z = 1 and Z — 2. Indeed, the same is true 
for higher Z, for Z between 3 and 8 the difference between 
C/ and Vefj for the Is atomic orbital, obtained from the 
appropriate ionization energies |l| is respectively , in eV, 
4.22, 4.22, 4.21, 4.19, 4.15, 4.05. We conclude that for Is 
orbitals one has approximately 



\J,fi{Z) = U{Z)-A.\eY. 
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for a large range of Z. A similar reduction from bare V to 
effective U will be found for other atomic orbitals. This 
reduction occurs because when the second electron is 



added to the singly occupied orbital the state of the two- 
electron system is not the doubly occupied single electron 
orbital. Rather, the orbital will 'expand' to reduce the 
Coulomb repulsion between electrons, and furthermore 
the two electrons will develop angular correlations. This 
effect is of course well known in atomic physics, and in its 
simplest form is approximately described by Slater's rules 
for the shielding constants 0: when another electron is 
added to an atom, the effective Z for the electrons in the 
same shell is reduced. This effect is however ignored in 
the ordinary Hubbard model 

A simple way to describe this physics is by introducing 
coupling to a fictitious local boson displacement coordi- 
nate qi for atom i that modulates the Hubbard U: 



U{qi) = U + aqi 



(4) 



that will relax when double occupancy occurs. As the 
simplest model we describe the boson dynamics by a har- 
monic oscillator of frequency ujo — {K/My/"^: 



H, = ^ + -Kqf 
2M 2 



(U + aqi)nii;nii 
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and the effective on-site repulsion is found by completing 
the squares as Uef f — U — a-/2K. The equilibrium 
position of the boson is gt = O for the orbital empty or 
singly occupied, and qi = —a/ K for the doubly occupied 
orbital. 

The reader may argue that the ordinary Holstein 
model ||] 



H^ = 
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+ n^li + OLqi{nii +ni{) + Uui^nii (6) 



will also describe a reduction of the bare [/ to a Uef / — 
U — a^/K. However, in contrast to Eq. (5), Eq. (6) also 
describes dressing of electrons in smgZy-occupied orbitals 
by the boson degree of freedom. That is not the physics 
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we are trying to describe here: without electron-phonon 
interactions or coupling to atomic electrons in other or- 
bitals, the electron in the singly occupied orbital should 
be undressed. 

In terms of boson crcation and annihilation operators, 
the Hamiltonian Eq. (5) is 

Hi = ujoala.1 + [U + gwo(aJ + ajjrii-fnii (7a) 



Ueff^U-wog'' (7b) 

with g = a/ {2KuJoy^^. The boson degree of freedom de- 
scribes the electronic excitation of an electron when a sec- 
ond electron is added to the orbital. Hence the frequency 
ujo is related to the excitation energies of the atom, and 
we expect 



H, 



l-b 



g{n)uJo{al + aj)(ni| + n^) 
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that is, an ordinary electron-boson coupling with a cou- 
pling constant that increases monotonically with band 
filling. Hence, as the usual electron-phonon interaction, 
it will give rise to an effective mass enhancement and 
a quasiparticle weight reduction, which increase as the 
band filling and hence the effective coupling g{n) in- 
creases. There are however other dynamical effects of 
Eq. (9) that are lost in the mean field treatment Eq. 
(10). 

To study Eq. (9) we perform a generalized Lang-Firsov 
transformation on the fermion and boson operators 



= Xi, 
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with c a constant of order eV, since the excitation en- 
ergies in an atom scale with the square of the nuclear 
charge. From Eqs. (3) and (7b) we conclude 
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with c' — A.leV/c. For a lattice system where an electron 
hops from site i to site j with hopping amplitude the 
Hamiltonian is thcn 

H = - ^ ti.j{c\„Cj„ + h.c.) -f ^[C/ + gujo{a\ + a,,)]ni|n,x 
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The physics described by Eq. (3), that led to the 
Hamiltonian Eq. (9), is a ubiquitous phenomenon, origi- 
nating in the fact that the spacing between atomic energy 
levels is always smaller than the strength of the Coulomb 
repulsion between electrons in a given orbital. Its quanti- 
tative importance is determined by the magnitude of the 
ionic charge Z, as discussed below. Other models, with 
auxiliary spin degrees of freedom |^ or with only elec- 
tronic degrees of freedom |^ , can be constructed contain- 
ing the same physics. The model Eq. (9) is a particular 
case of the generalized Holstein models discussed in ref. 

By construction, the model Eq. (9) is not electron-hole 
symmetric, even if the band structure defined by t.y is. In 
particular, a single electron in the empty band does not 
interact with the boson ficld at all. For a few electrons in 
an empty band the effect of the bosons is negligible if the 
bare Coulomb repulsion U is appreciable, as the proba- 
bility of double occupancy will be small. In contrast, a 
single hole in the fuU band interacts most strongly with 
the boson field. Treating the four- fermion term in mean 
field, the electron-boson part of the Hamiltonian Eq. (9) 
is 



0.1 - grii^nii 



(Hb) 



and the Hamiltonian Eq. (9) becomes 



(12) 



with Uef f given by Eq. (7b) . The ground state expecta- 
tion value of the operator is 

< X,, >o= e-(0'/2)».- = 1 + (^ - l)n,,_„ (13a) 



S = e-^'/2 



(13b) 



The part of the fermion operator Eq. (lla) associated 
with ground state to ground state transitions of the bo- 
son field is the coherent part of the operator, the quasi- 
particle. We have then 

C, = |0 X 0| [1 + (5 - l)n.,-.]5.. -f C™'' (14) 

The incoherent part of the operator 



c™™'' = I n. 
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describes processes where the boson field makes transi- 
tions to and from excited states \l >,l ^ O, which only 
take place if hi^-^ = 1, that is if the orbital is occupied 
by another electron of opposite spin. 

The quasiparticle weight in this model is, from Eq. 
(14) 



z{n) = [l + -{S-l) 



(16) 
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and decreases monotonically with electronic band filling 
n, O < n < 2, so that quasiparticles become increasingly 
dressed as the band filling increases. The factor S is the 
overlap niatrix element of the oscillator ground states 
with and without site double occupancy [Q| , and gives 
the quasiparticle weight for a hole in the filled band (n = 
2 in Eq. (16)). According to Eq. (8b), as the ionic 
charge Z decreases S' decreases rapidly, implying that 
hole quasiparticles become increasingly incoherent. 

We can estimate S from first principles for a hydrogen- 
like ion. In the Hartree approximation, S will be given by 
the overlap matrix element of the electron wave function 
in the presence and in the absence of another electron in 
the orbital: 



S 



< ipis\tpis > 



(I 5_^3/2 

m Z I 
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with Cpig the Is orbital Eq. (Ib) with Z replaced by 
Z = Z — 5/16, as appropriate for the Hartree wavefunc- 
tion. If we use the more accurate Eckart wave function 
P.pii[ ^Eck for the two-electron ion, which incorporates 
radial correlations, we can estimate S from the square 
root of the overlap matrix element of the Eckhart wave 
function and the wavefunction of the two electrons in the 
Is orbital, as 



S' = a/I < (pis'Pls\i^Eck > 



2f{Z, Zi)f{Z,Z2 
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f (Z Z') 



f{z,,Z2)y/^ 

(18a) 
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Here, Zi and Z2 are the two orbital exponents for the 
Eckart wave function obtained by minimization of the en- 
ergy, which are found to be approximately Zi — 1.14Z — 
0.105, Z2 = 0.905Z - 0.622. The Eckart wave function 
becomes unstable for Z < 0.93. 

Figure 1 shows the dependence of S on the ionic charge 
Z. We expect the Eckart wavefunction to underestimate 
and the Hartree wavefunction to overestimate the value 
of S for a given Z. Figure 1 also shows S for the dynamic 
Hubbard model, Eqs. (13b) and (8b), for a value of c' 
that matches the Eckart wavefunction results for large 
Z. While the three curves are different, the qualitative 
behavior is the same, showing a decrease of S as the 
ionic charge Z decreases, which describes the increased 
incoherence of single hole carriers in the system as the 
ionic charge decreases. 

Upon replacement of the operator form Eq. (14) in 
the Hamiltonian Eq. (12), and ignoring the incoherent 
part of the operators, we obtain the effective Hamiltonian 
describing propagation of quasiparticles: 

Hqp^~Yl *y (^Lcja + h.c.) + UeffHi^nii (19a) 



where the hopping amplitude now will depend on the oc- 
cupation of the two sites involved in the hopping process: 

tlj = ^ufl + (S - l)(ni.-cr + «j,-(t) + (S - l)^nj,_^nj,_cr] 

(19b) 

The effective quasiparticle Hamiltonian Eq. (19) will ac- 
curately describe the low energy physics of the fuU Hamil- 
tonian Eq. (12) in the strong coupling regime, where 
the Lang-Firsov approximation to the Holstein model be- 
comes accurate. Even in that regime however the high 
energy degrees of freedom described by the fuU Hamil- 
tonian Eq. (12) still play an important role, in ensuring 
that various sum rules that are violated by the low energy 
Hamiltonian Eq. (19a) are satisfied. 

The hopping amplitude for an electron of spin a when 
there are electrons of opposite spin at both sites i and j, 
i. e. Ui^-a- — fij.-a = 1 is, from Eq. (19b), t2 — S'^tij, 
and when there is onc other electron at the two sites, i. e. 
ni,-a + ^j,-cr = 1, ^1 = Stij. If S « 1 then t2 « ti. 
In the limit where t2 can be neglected with respect to ti 
it becomes obvious that the Hamiltonian Eq. (19) leads 
to pairing of two holes in a fuU band: a state where the 
two holes are separate describes essentially localized holes 
with zero energy, since f2 ~ 0; a lower energy state results 
from a linear combination of states where the two holes 
are on the same or nearest neighbor sites, with energy 
of order Cpair ^ — 2p x t\/Ueff with p the number of 
nearest neighbors to a site. More generally, an exact 
criterion can be found for the parameters that will yield 
a bound state for two holes in a filled band described by 
the Hamiltonian Eq. (19) [O: 



U, 



D 



< 1-S^ 



(20) 



where D is the bare bandwidth of the band defined by tij . 
The same criterion is found for existence of superconduc- 
tivity in the dilute carrier concentration regime in BCS 
theory , a regime where BCS theory is expected to be 
accurate, and quantitatively close results are obtained 
from exact diagonalization of small systems (l3| . As the 
atomic charge Z decreases, the coupling g increases (Eq. 
(8b)), hence S decreases (Eq. (13b)), and furthermore 
Ueff decreases (Eqs. (3) and (la)). Both of these ef- 
fects are in the right direction to satisfy the condition 
Eq. (20). Furthermore, if the interatomic distance is de- 
creased, and the bare bandwidth D will increase, again 
in the direction of satisfying the inequality Eq. (20) . We 
conclude that superconductivity induced by this physics 
will occur in systems where the ionic charge Z is small, 
the interatomic distances are small, and the Fermi level 
is close to the top of a band. The parameter regime 
where the Hamiltonian Eq. (9) yields superconductivity 
has not yet been established, however numerical studies 
of the similar spin- 1/2 Hamiltonian Q have shown that 
pairing survives well beyond the parameter regime where 
the mapping to the effective Hamiltonian Eq. (19) is ac- 
curate. 
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The dressing of hole quasiparticles decreases in this 
model as the hole concentration increases, and corre- 
spondingly the hole quasiparticle weight increases, as 
seen from Eq. (16) or its equivalent in hole represen- 
tation 



(21) 



with Uh — 2 — n, T — I/ S — 1. In the normal state, this 
should be seen in the one-particle Green's funetion as a 
transfer of spectral weight from the incoherent part, de- 
scribing high energy excitations of the boson field on the 
scale of ujo, to the quasiparticle peak as the hole dop- 
ing increases. Correspondingly, a transfer of spectral 
weight in optical absorption should be seen, from high 
frequency absorption on the scale of loq to intra-band 
Drude absorption with plasma frequency determined by 
t2{nh) = ^2 + nh{ti — t2). From the operator relationship 
for hole quasiparticles derived from Eq. (14) 



(22) 



it is seen that the expectation value < dla'^fa > 
quires contributions from anomalous expectation values 
< c^„c% >. This will cause transfer of spectral weight 
from high to low frequencies in the one- and two-particle 
Green's functions for fixcd hole doping when the system 
goes superconducting, which should also be observable in 
photoemission and optical absorption experiments . 

Non-perturbative theoretical techniques such as dy- 
namical mean field theory jl4j and DMRG [|l5| should 
be able to establish the parameter regime in the Hamil- 
tonian Eq. (9) where the physics of superconductivity 
through undressing described above takes place. First- 
principles quantum chemical and density-functional cal- 
culations should be able to relate the parameters in Eq. 
(9) to real materials. The physics discussed here predicts 
that superconductivity is favored in systems where con- 
duction is through holes in nearly fiUed bands, through 
conducting structures that are negatively charged (small 
ionic charge Z), and that such systems will show undress- 
ing of hole carriers both when the system is doped with 
holes and when it goes superconducting. These effects 
will be most apparent when Z is small and T is large, 
which also leads to high Tc- The facts that most su- 
perconductors show hole carrier transport in the normal 
state , that both MgB2 and the high Tc cuprates have 
holes conducting in highly negatively charged substruc- 
tures {B^ and {Cu02)^ planes), and that the high Tc 
cuprates show evidence of 'undressing' upon hole doping 
and upon going superconducting [^-^ T), suggest that 
the physics of the dynamical Hubbard model may have 
something to do with the physics of superconductivity in 
real materials. 

In the regime most favorable for superconductivity, i. e. 
hole conduction in a system with large T, the model 
Eq. (9) predicts that incoherent excitations on elec- 
tronic energy scales (wg) will be induced when holes hop. 



Hence one would expect non-thermal high frequency ra- 
diation to be generated when a dc current circulates in 
the normal state. This electroluminescence should 
be most pronounced and easily observable at the posi- 
tive counterclectrode-sample boundary, where holes are 
injcctcd into the sample. The intensity of the radiation 
should correlate with the magnitude of the dressing and 
hence of Tc, and the frequency distribution will give In- 
formation on the scale of electronic excitation energies 
involved in the undressing process that leads to super- 
conductivity. 
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FIG. 1. Parameter S versus ionic charge Z in the Hartree 
and Eckart approximations, and in the electron-boson model 
with c' = 0.2 (Eqs. (8b), (13b)). The deviation of S from 
1 indicates the degree of quasiparticle dressing for an almost 
fuU band. 
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